Casimir Force between Vortex Matter in Anisotropic and Layered Superconductors 
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We present a new approach to calculate the attractive long range vortex-vortex interaction of 
the van der Waals type present in anisotropic and layered superconductors. The mapping of the 
statistical mechanics of vortex lines onto the imaginary time quantum mechanics of two dimensional 
charged bosons allows us to define a 2D Casimir problem: Two half-spaces of (dilute) vortex matter 
separated by a gap of width R are mapped to two dielectric half-planes of charged bosons interacting 
via a massive gauge field. We determine the attractive Casimir force between the two half-planes 
and show, that it agrees with the pairwise summation of the van der Waals force between vortices 
previously found by Blatter and Geshkenbein [Phys. Rev. Lett. 77, 4958 (1996)]. 

PACS numbers: 74.60.-w, 74.72.-h, 74.90. +n, 74.25.Dw 



I. INTRODUCTION 



Within the Shubnikov phase of type II superconductors 
the applied magnetic field enters the sample in the form 
of flux lines. The standard mean- field type calculationtl 
shows that in an isotropic material two straight vortices 
repel at all distance scales, with an interaction strength 
T^rcp = '2eoKo{R/X), Eq = (<I>o/47rA)^ being the basic 
energy scale in the vortex matter, Kq is the zero-order 
modified Bessel function, R is the inter-vortex distance, 
and A the magnetic penetration length ($0 = hc/2e de- 
notes the flux quantum). However, it has recently been 
showna that in layered and strongly anisotropic super- 
conductors the thermal fluctuations of the flux lines give 
rise to a long range attraction Kdw ^ ~{T/d){X/R)^ of 
the van der Waals type between the vortices, where T de- 
notes the temperature and d is the interlayer separation. 
The strongly fluctuating and layered high temperature 
superconductors are particularly well suited to exhibit 
this attractive component in the vortex-vortex interac- 
tion. Alternatively, the attraction is induced through 
static vortex distortions due to an underlying pinning 
landscape, an effect recently studied by Mukherji and 
NattermannB and by Volmer et ala 

Following a suggestion of Nelsoni, the statistical me- 
chanics of vortices can be mapped to the imaginary time 
quantum mechanics of two-dimensional (2D) bosons. 
The particular interaction between the flux lines renders 
the bosons charged (with a charge screened on the scale 
of the London penetration depth A). This type of long 
range interaction can-be formulated in terms of a massive 
gauge held theoryQ^H. Within the resulting 2D massive 
electrodynamics, the vortex matter acts as a dielectric 
medium and we can define a Casimir problem, see Fig. 1: 
Under the vortex boson mapping two half-spaces of 
vortices separated by a gap of width R act as two dielec- 
tric planar media which attract each other via a Casimir 
force. 




FIC. 1. Ceometry used in the calculation of the Casimir 
force between two half spaces of vortex matter separated by a 
gap of width R. Shown are two parallel cubes of size x 
at a distance R with vortex density p (Lz — > 00). We choose a 
coordinate system with the z-axis parallel to the vortex lines 
and the a;-axis perpendicular to the planes of the gap (dark 
shaded planes). The statistical mechanics of the flux lines 
is mapped to the imaginary time quantum mechanics of two 
dimensional (2D) bosons (solid circles). The cut through the 
vortex lines (lightly shaded planes) defines a time slice of the 
2D bosons. The Casimir energy produced by the dielectric 
action of the two media is calculated by summing over the 
ground state energies of the fluctuating massive electromag- 
netic field mediating the interaction between the 2D bosons. 
Using the technique of pairwise summation (dashed line) valid 
for rarefied media (i.e., p 0), we can relate the macroscopic 
Casimir effect with the microscopic van der Waals attraction 
between the particles involved. Picture not to scale. 

In the present work, we determine the Casimir force 
between two dielectric half planes of charged bosons. For 
dilute media, the macroscopic Casimir force can be re- 
lated to the microscopic van der Waals force between the 
media's constituents via pairwise summation. Here, we 
present a derivation of the van der Waals force between 
vortex lines via this alternative route, calculating first 
the Casimir force between two dielectric planar media in 
the boson picture and then reconstructing the van der 
Waals force in the reductionist way. 

Since the discovery of the Casimir effecJ^ in 1948, sev- 
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eral hundred papers have dealt with this phenomenon, 
dissem«ia,ti«g its fascination into many branches of 
physicgl3il3. Casimir forces between macroscopic bod- 
ies are a quantum effect caused by a shift in the zero- 
point energy of gauge field fluctuations such as the elec- 
tromagnetic one. The Casimir effect is bound to a num- 
ber of system properties, such as topology and dielectric 
permittivity, and the reduction to an analogous van der 
Waals interaction is not always possible. However, the 
interpretation in terms of a van der Waals attraction is 
possible for the case of rarefied media and appropriate 
geometries, such as the parallel plate setuptZl. In our 
derivation of the van der Waals force from the Casimir 
effect we will make use of these special conditions. 

In the following, we discuss the relationship between 
the Casimir effect and the van der Waals interac- 
tion within a path integral formulation (Sec. II) be- 
fore deriving the appropriate action for the 2D charged 
bosons from the London functional describing the vor- 
tices (Sec. III). In section IV, we briefly review the deriva- 
tion of the van der Waals force in the original vortex lan- 
guage and then proceed with the calculation of the 2D 
Casimir effect in section V, the main section of the paper 
containing the new results. 



II. CASIMIR VERSUS VAN DER WAALS 

We consider two parallel material slabs made from fluc- 
tuating dipoles and separated by a vacuum gap. Sum- 
ming pairwise over all microscopic van der Waals in- 
teractions between the dipoles provides the macroscopic 
Casimir interaction between the slabs. On the other 
hand, we can determine the dielectric properties of the 
individual slabs as produced by the fluctuating dipoles. 
The specific boundary conditions due to the dielectric 
properties of the slabs influence the spectrum of the elec- 
tromagnetic field confined in between. The change in 
the spectrum as a function of the separation of the slabs 
produces the Casimir force. This analogy between the 
Casimir and the van der Waals force is transparently 
brought out within a path integral formulation, see Fig. 2: 
Assume the system under consideration can be described 
by an action 5[a,j] depending on the gauge field a and 
a particle current j. Carrying out the partial integration 
in the partition function Z = J V[a]V[i] e^'^U''^!/'^ over 
the matter field j or over the gauge field a, we obtain an 
effective action SoS describing the conjugate field alone: 
From the effective action SoS [a] describing the gauge field 
we can derive the Casimir effect, while the interaction of 
the particle currents as described by 5off[j] will give us 
the van der Waals attraction. 

The Casimir- and van der Waals forces then can be 
related to one anotjher via the pairwise summation of the 
interparticle forcestj: Consider two d-dimensional homo- 
geneous macroscopic dielectric bodies of density p with 
parallel interfaces separated by a distance R. The two 



bodies attract one another due to a microscopic particle- 
particle interaction V^dw = Ar~", a > d, of the van der 
Waals type. The interaction energy can be written in the 
form 

pOC POO f' 

U{R)=L''-'p'' dx dx' d^-iyKdw(r), 

JR Jx JE 

where E' is a hypercube of size L'^"^ parallel to the in- 
terface and r = y^r^ + j/^, with y the — I dimensional 
in-plane coordinate, while x is the coordinate along the 
direction perpendiculax to the plane. We then find for the 
Casimir force densityll3 / —dRU{R)/L'^^^ the result 



f = P^J^ / d'^-iyKdw(r) 



^ .„2_(.-i)/2 r[(l + a-rf)/2] 1 1_ 

T{a/2) a-dR'^-" 



(1) 



The result (^ then allows to infer the parameters A and 
a, characterizing the van der Waals interaction Vvdw, 
from the macroscopic Casimir force density /. 
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FIG. 2. Schematic relation between the Casimir effect and 
the van der Waals attraction. 



After mapping the thermally fluctuating vortex matter 
to a system of 2D quantum charged bosons, we will ar- 
rive at an action A, a] with two gauge fields a and A, 
see Eq. (|^). The Casimir effect we are interested in here 
is the one produced by the fake gauge field a: The inte- 
gration over the matter field j will produce the dielectric 
properties of the media, while the integration over the 
physical gauge field A renders the fake field massive and 
hence will always exponentially confine the Casimir force 
to finite distances. 



III. BOSON REPRESENTATION OF VORTICES 



We start from the London free energy in an isotropic 
superconductorEI 



.F[B] = i- j dV[B2 + A2(VxB)2], 



(2) 
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with B = V X A the magnetic field and A denoting the 
London penetration depth. In order to acquaint for the 
vortices, we add the current term — (j • B)<i>o/87r with 



j = (J,p)=^t,,(z)<5(2)(R_R^^(^)), 



(3) 



where r = (R, z), R = {x,y) E K^, the coordinates 
R^j(z) denote the position, the vectors = (dzH^iz), 1) 
the direction of the vortex lines, and $o — hc/2e is the 
unit of flux. Ignoring screening, the interaction between 
the vortex lines is long ranged and thus can conveniently 
be expressed through a mediating gauge field a. As usual, 
we introduce the gauge field a as an auxiliary field such 
that /X»aexp{-/3J^'[a])p= exp{-f3J^), with f3 = 1/T 
the inverse temperatures (we set the Boltzmann con- 
stant fcs to unity and fix the gauge through the condition 
V-a = 0), 



.f'[a,A,j] =Jd' 



(4) 



— (Vx A)2 + — (Vxa) 



25^ 



where g"^ = 47reo and the line energy sq ~ ($o/47rA)^ is 
the basic energy scale in the problem. In (|^), we have ac- 
counted for screening by introducing back into the model 
the real gauge field A. 

The anisotropy of uniaxial layered materials is most 
conveniently introduced through, a- i ra scaling of the scalar 
fake magnetic field 6 = (V x a)^^^, 



(V X a). 
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(V X a),. 



where the anisotropy factor e — ^jmablTnc is deter- 
mined through the effective masses perpendicular (mc) 
and parallel (mab) to the a6-plane (the fake electric field 
e_L = (cj,,— Ba;) = —(V X a)xy remains unchanged; the 
subscript 'xy' identifies the planar component of a vec- 
tor, a = (a^ryjOz)). 

We map the statistical mechanics of the vortex system 
to the imaginary time quantum mechanics of 2D bosons 
through the replacementgatSS z ^ t (imaginary boson 
time), hP (the boson's Planck constant), T' ^ S 

(the boson action), and 1/e ^ c (the light velocity in the 
boson system). The boson partition function is given by 



Z = y I?[{R^}]X>[A]2:»[a]e-'5[^^^>^^^'^l/'"''', 
with 5 = iSo + 5int and 
5o[{R4]= /'dr^[^(9.R4r))^-/. 

5i„t[{R^},A,a] = ^AtA^R 



(5) 



(6) 



*a- lJ-^(Vx A) ) (7) 
^(VxA)2 + ^((Vxa)^^ + l(Vxa)? 



Here, yiP = H^o/iir—eo ln(A/^) is the chemical potential 
of the bosons, p = B/^q is the boson density, and m — Eq 
is the boson mass (we have introduced a term — BH/47r 
to account for the external magnetic field H producing 
the vortices in the superconductor; ^ denotes the planar 
coherence length). In the thermodynamic limit, we have 
Lz — > oo, corresponding to = HP /Lz — > 0, i.e., we 
are interested in the ground state physics of the boson 
system. Note that the free boson action iSq contains the 
bare mass m due to the vortex core energy. The retarded 
self-interaction of the bosons via their gauge fields then 



produces the mass renormalization m 



siikz), 



where ei(kz) is the dispersive line tensionEl of the vortex 
lines. 



IV. VAN DER WAALS INTERACTION 



To set the stage, we briefly review the derivation of 
the van der Waals interaction as presented in Ref. H; this 
will allow us to fix some flaws in the previous derivation 
and will provide us with a check on the results for the 
Casimir force derived later. 

For a simple qualitative analysis we consider two vor- 
tices in a layered superconductor (with layers separated 
by the distance d). Ignoring the coupling between the 
layers, the fluctuating pancake vortices interact via a 
logarithmic potential V{R) « 2eoln(i?/A). Second or- 
der perturbation theory then provides us with a van der 
Waals interaction — Kdw ^ {T/d){X/Ry, the energy 
scale T/d being set by the driving action of the ther- 
mal fluctuations. At long distances R > d/e a finite 
interlayer coupling changes this result as the Josephson 
interaction reduces the cutoff 1 /d on the kz modes to the 
new value 1/eR; the van der Waals interaction crosses 
over to Vvdw ^ —{T/eX){X/R)^^. The two results have 
their analogue in the van der Waals attraction between 
neutral atoms, where the interaction potential exhibits a 
crossover-from a at short to a behavior at large 
distancesE3. 

We proceed with the derivation of the long range van 
der Waals interaction between vortex lines/bosons. Fol- 
lowing the scheme in Fig. 2, we integrate over the gauge 
fields A and a in the partition function (|5|) and obtain 
the effective current-current interaction 
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d3rdVj„(r)F^f(r-r')jMr) 



(8) 



(we return to the more natural statistical mechanics no- 
tation in this section. Sen ^ J^)- Inserting the expression 
® fw the currents, we can cast Eq. (||) into the standard 
forriM 



£o \ ^ 
2 ^ 



(9) 
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with the vortex positions = (Rp,z) and the inter- 
action potential V^, conveniently expressed within a 
Fourier representation, 



(2^ 



^ik(r— r') 



K";(k), (10) 



with 



1 + A2fc2 



Sa0 — 



(11) 



Here k — (K, fcz) — {kx,ky,kz), K_l = {ky,—kx) and 
\c — X/e, S,c — The free energy (^) can be split into 
the self energy part {p — v) and the interaction part 
(m 7^ Restricting ourselves to two vortices a 
distance R apart, we obtain 



dzidz2 e^K[R.+ui-u2]g»fe.(^i-^2) 



4^7 (27r)3' 
X [Kf(k)+ii.(zi)i2/3(^2)K^'(k)], 



(12) 



where we have split the vortex positions {Ri,R2} into 
a mean field part {R?,R^} = {0,R} and a fluctuating 
part {ui,U2}. Up to a constant, the free energy is given 
by 

F{R) = -T\nZ{R) = -rin(exp[-/3JC-'"*])o. (13) 

The average (. . .)o has to be taken with respect to the 
self-energy J-'^ of the free vortices. Performing a cumu- 
lant expansion, we obtain the effective vortex-vortex in- 
teraction in the form L.V^s « {J''"^^)o - [{{^"'^?)o - 
(JF"^*)q]/2T, where denotes the sample thickness. 
Splitting into longitudinal (to the induction, JFn) and 
transverse parts, the longitudinal interaction pr 

duces the standard repulsive vortex-vortex intcractior 
to lowest order in u^,. 



I 



Kcp(i?) - 2eoKoiR/X), 



(14) 



while higher orders in merely renormalize the pre- 
factor. The. transverse part produces the van der Waals 
interactiono 



dW 



(15) 



Using the decomposition 

{tla{~kz)t2p{kz)tia' {-k'^)t2t3' ik'J)o 
= 2TTLz5{kz - k'JSaa'SpfS'- — ^"T^, 



we can reduce the average in (|15|) to the simpler form 
dfc. 



647r2 ./ 27r 



[V^1}{R,k,)Y{t'{k,))i (16) 



with the partial Fourier transform 



v:'^{R,k. 



{2n 



-K^'(K,fc,)e 



iKR 



(17) 



In strongly anisotropic and layered material the single 
vortex mean squared amplitude of fluctuations (t'^{kz))Q 
is limited by the electromagnetic interaction through the 
dispersive elasticity ei{kz), 



2T 

£i{kz) ' 



with 



ei{kz) 



£0 



2A2fc; 



■In(l + A2fc2). 



(18) 



(19) 



The evaluation of the partial Fourier transform (|r^) is 
carried out in Appendix A and making use of the results 
(Al) and (A2) in the limit e = 0, we find the van der 
Waals interaction in the decoupled limit 



4£o 



T fx 



vdW 



ln2(7rA/d) dea \R 



(20) 



The interlayer distance d provides the large kz cutoff on 
the thermal fluctuations of the individual lines. 

In the continuous anisotropic case the single vortex 
elasticity is still dominated by its electromagnetic con- 
tribution. At intermediate distances A < R < d/e we 
recover again the result (po| ) as the Josephson coupling 
is not effective yet. At large distances A, d/e < R < 
Ac, however, the interlayer coupling becomes important 
through cutting off the kz modes at l/eR and we find 
the result 



•'vdW 



(17l7r/256)eo T /A 



(21) 



These results coincide with the ones found by Blatter 
and Geshkenbeino up to the numerical prefactor as they 
missed a factor 2 in Eq. (p^), as well as a term in the 
partial Fourier transform which contributes to the result 
( pT| ) to order 0{s~^). The same mistakes show up in 
their results for the attraction of the vortex to the sam- 
ple surface, Eqs. (19) and (20) of Ref. |. We give the 
correct expressions here: 



vdW 



£o/2 T fx 
ln(7rA/d) dea \R 



for the decoupled case, and 
(3/16)£o 



v: 



T 



vdW 



ln(7rA/ei?) eAeo \R 



X 



(22) 



(23) 



in the continuous anisotropic situation. 
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V. CASIMIR FORCE 

In this section, we derive the Casimir force between 
two parallel cubes of vortex matter and then show, via 
the method of pairwise summation, that the results co- 
incide with those found by means of the van der Waals 
approach. 

Again, we begin with a simple analysis based on dimen- 
sional estimates, instructing us what to expect. Consider 
two d-dimensional hypercubes of size L"* separated by a 
distance R ^ L, see Fig. 1. We wish to determine the 
Casimir force / per unit area acting between the inter- 
faces. Being due to the zero point fluctuations of the 
gauge field, the only relevant dimensional quantities en- 
tering the expression for the force are Planck's constant 
fi,^, the velocity of light c, and the distance R between 
the slabs. The combination h^c/R'^'^^ then provides us 
with a dimensionally correct expression for the force; in 
two dimensions we then find / ~ h'^c/R^ for the re- 
tarded Casimir effect. In the non-retarded case, c ^ oo 
and we have to replace the combination c/ R by the fre- 
quency scale w where the matter becomes transparent 
(for the vortex matter in a layered superconductor this 
'frequency' is given by the layer separation d, lu ^/d). 
Dimensional estimates then give us for the force per unit 
length in two dimensions / ~ h'^uj/R^. In order to de- 
rive the corresponding microscopic inter-particle poten- 
tial from these estimates, we apply the method of pair- 
wise summation and obtain the result V^dw 1/sR^ in 
the retarded case and V^dw oc in the non-retarded 

limit. These simple estimates produce the correct power 
laws for the van der Waals potential; in order to find the 
correct sign and the complete prefactors, we have to go 
through the detailed analysis below. 



A. Formalism 

The Casimir energy is given through the differ- 
ence in the sum rOwec all cavity modes minus the free 
field contributiorO'CJ. While the calculation is rather 
straightforward for the simple case of a metallic cavity, 
an appropriate formalism has to be set up to treat more 
general configurations involving dielectric media. 

We start from a quadratic Lagrangian density C = 
aG~^a in d dimensions, with a(x, t) a vector boson field 
with r components and G the Green function matrix. 
The partition function Z is expressed through the usual 
imaginary time path integral formalism 



Z— I?[a] exp 



'X I dr C[a\ 
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[det G] 



and the free energy F — —T\nZ is given by 



F = -TrlnG-^ 
2 



^1/2 



(24) 



(25) 



We consider the classic parallel-plate geometry with 
isotropic media separated by a gap of width R along the 
Xd-direction. For each polarization j/ e {1, . . . , r} we de- 
termine the eigenstates obeying the boundary conditions 
at the {d — l)-dimensional hypersurfaces placed at and 
R. The individual modes are characterized by their po- 
larization h', the (Matsubara) frequency i^s = 2TTsT/h, 
s e Z, the transverse wave vector q e R''^^, and the lon- 
gitudinal wave vector A;„, n e N. Given the set {v, ^s, q), 
the discrete wave vectors fc„ are obtained as the solutions 
of the boundary condition D^g^y{k) = 0, where {a} is a 
short hand for the other indices {i^,£,s,ci} (the function 
D^^j(k) derives from the determinant associated with 
the set of boundary conditions; for the classic Casimir 
effect, D{k) ~ sin(fci?) with R the distance between the 
metal plates). Expressing the trace through all the in- 
dices we have 



F : 



TL 



■EE' 

u=l s=Q 



^ InG^Jj(fc), (26) 
{c}(fc)=o 



where the prime on the sum indicates that we count the 
s — term with a weight 1/2. In ( p6| ) we have made 
explicit the boundary condition at the hypersurfaces sep- 
arating the media from the gap. 



Imk 



X 



Rek 



FIG. 3. Contour and distribution of poles in the complex 
fc-plane as they appear in the expression (|27|). The poles 
produced by the boundary condition D{k) = are marked 
by the solid dots on the real axis. The zeros of g{k) are de- 
noted by the crosses. All zeroes are distributed symmetrically 
with respect to the origin, a consequence of the symmetry 
G{k) = G{-k). 

We make use of Cauchy's theorem to perform the sum 
over the longitudinal momenta fc„: We first rewrite the 
sum in the form 

E ln[G-+y]=/d, J: 

Next, we let the boundary condition D^a}{k) — pro- 
duce the desired sum in a Cauchy loop integral in the 
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complex k plane, see Fig. 3. With g{k) = (G|Jj(fc) + y) 
we have 



1 



dkPjk) 1 
D{k) g{k) 



(27) 



y — 

^ g(k) 



E 

g(fe)=0 



D{k) dk9{ky 



where it is understood that both D and g further depend 
on the other indices {a} as well as on the parameter y in 
the case of g. We define the function 'w{{a]\ y) as the ze- 
roes of the expression g{a}ik) = G^^y[w{{a}; y)]+y ^ 0. 
Taking the derivative of the last equation with respect 
to y, dyg = dyG^^ + 1 = dkG^^dyW + 1, and using the 
result in (127^), we obtain the desired result 



1 - = 2|- In [«;({«}; y)], (28) 
y oy 



where we have used that the Green function exhibits the 
proper asymptotics G{k) ^ 0{k^) and the symmetry 
G(fc) = G{—k). Inserting the sum on k, Eq. (gsj) with 
y — Q, back into the free energy expression (|2^) we obtain 

:F = TY,\nD{^}[w{{a})], (29) 



where 



denotes the sum/integral over all modes {a} = [v^ £,s,^)■ 
The above formal manipulations require a regulariza- 
tion such that D[^s — ^ oo] 1. We carry out the 
appropriate renormalization by the subtraction of the 
free energy of the free field without boundary condi- 
tions, corresponding to replacing the partition function 
through the ratio Z — > Z/Zq, with Zq the partition 
function of the free field. For the latter, we define the 
'boundary condition' Dq in such a way as to produce 
the fraction R/L oi the free field energy in the volume 
L'^ via the expression (pof); for the classic Casimir prob- 
lem, Do(fc ± i5) cx limL^oo{sin[(fc ± i5)L\}^/^. With 
£)rcn _ D/Dq we finally obtain the regularized free en- 
ergy of the Casimir problemO 



({«})]. 



(30) 



{a} 



We illustrate the use of this formalism with a brief 
derivation of the classic Casimir result, the attraction 
between two parallel metallic plates of size L x L and 
a distance R apart. The dispersion relation for the 
free electromagnetic field is G^^{^s,<l,k) = + 
+ k"^, such that the function w\y=o takes the form 

w(^s,q) — iy/S.s/c'^ + q^- The boundary condition re- 
quiring the fields to vanish at and R can be cast into 



the form D{k) = sin(fci?) = 0, producing the modes 
with a longitudinal k vector fc„ = 27rn/i?. The 'bound- 
ary condition' describing the free field takes the form 
Dn{k±iS) = (=Fl/2i)exp(=F«fci?), such that £''■•="(^±^(5) = 
1 — exp(^2ifci?), where S 0+. Combining the results 
for D^""^ and w, we have to carry out the mode summa- 
tion over the logarithm of 



i?'™[u'(6,q)] = 



1 - cxp -2R 



£2 



n 2 



(31) 



where taking the square accounts for the two polarization 
modes of the electromagnetic field (the single polariza- 
tion mode at A: = has been properly taken into account 
here). The force density / between the plates is given 
through the derivative of the energy (|30| ) with respect to 
R, 



I 



1 aj^'" 

l2 dR 



h 

ttR jo 



dqq 



2tt 



2RV¥/c 



exp(2i?VWcM^) - 1^ 



(32) 



where we have replaced the Matsubara sum through an 
integral over complex frequencies £. The remaining in- 
tegrals are easily evaluated by changing to polar coordi- 
nates q = p cos 0, £/ C|-p p sin </> and we obtain the classic 
result due to CasimirO 



/ 



TT^ he 

MOR^' 



(33) 



B. Application to 2D bosons 

In the following, we apply the above formalism to the 
2D charged bosons described by the action S — So + Sint 
as given in Eqs. (^) and (0). Following the scheme out- 
lined in Fig. 2, we first have to integrate over the cur- 
rents J ^ dr^i- We split the boson positions R/i(T) into 
the mean field part R°(t) and a fluctuating part u^(t), 
R^(t) = R°(t) -I- u^(t). In carrying out the integration 
over the boson positions R;i(T) we can ignore the fluc- 
tuating part u^(t) in the gauge field a as the latter is 
smooth on the scale of the amplitude u^(t), 

|pK]exp|-^|dr5:[^(9.u,)^ 



-*(a.u^,l)-a(R",r) 



(X exp 
where 



J d'RdT Qn(R)a2^(R) + *p(R)a,(R) 
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n(R) 



p(R) 



m 



and p(R)= V<5(2)(R-R") 



denote the polarizability and the density of the 2D 
bosons, respectively. Second, we have to integrate over 
the true gauge field A (we have fixed the gauge to 
V • a = 0). The resulting term (l/2(7^A^)a^ in the ac- 
tion then renders the fake gauge field a massive (physi- 
cally, this mass term expresses the finite range A of the 
interaction between the charged bosons/ vortices). Third, 
we introduce the free Green function matrix Go for the 
gauge field, in Fourier representation. 



Sxy 






,-1 



with {gxv)Jj = k'^Sa 

9r 



V/afj ^ "a/3 



(c2 - l)if2p^^(K), 



and the transverse projector Pap — 5ap — KaKp/K^. 
Combining the results of the above three steps, we arrive 
at the desired effective action for the fake gauge field a, 



5cfF[a] 



1 



Ixy 



"2.9^ 



9t 



1 



A2 



+ .g^n I a^ 

I- ipa^ 



(34) 



The Euler-Lagrange equations of this action determine 
the (imaginary time) field equations for the transverse 
and longitudinal parts of a (we remind that Ac = cA = 

A/£), 



9'p 



1 

a! 

^xy(yxy ■ i^xy) = (35) 



and 



xy 



1 

A2 



(36) 



with V xy = {dx, dy ). In the last equation of ( p5[ ) we have 
made use of the gauge condition V-a = \/ xy-axy+drCiT = 
and have ignored the time dependence in the longitu- 
dinal field component Ur, see below. 

The longitudinal field is generated by the source 
g^p/iiT = Eop and has no dynamics on its own. Placing 
two bosons at a distance R, a simple integration gives 
the vortex- vortex interaction (c.f., Eq. (14)) 



K-cp(i?) - 2eoKoiR/X). 



(37) 



Here, we ignore the fluctuations of the vortex position 
in the calculation of a,- as they merely produce a small 
renormalization of the prefactor in the repulsive poten- 
tial. Similarly, we can neglect such corrections in the 



calculation of the transverse modes below and set the 
r. h. s. of ( ^ ) equal to zero. 

The dynamical transverse field a.xy generates the 
Casimir force. Here, we have in mind a geometry as 
sketched in Fig. 1, with two parallel cubes of vortex mat- 
ter of density p, separated by a vortex free region of thick- 
ness R. We choose the x-axis to lie perpendicular to the 
'plates', the z-axis is directed along the vortices. Within 
the boson language we deal with two parallel planes of 
2D bosons lying in the xy-plane, the direction along the 
vortices now transforming to the imaginary time coor- 
dinate (below, we consider the limit — > 0, implying 
vortex lines of infinite length, Lz —^ oo). Note that in the 
present 2D case we have to consider only one polarization 
mode for the transverse gauge field Sixy, a consequence of 
the gauge condition V • a = 0. 

Going over to Fourier space we can cast (|3^) into the 
form (we remind the reader that c — q and k are the 
wave vectors along and perpendicular to the 'plates') 



e 



1 



1 

a! 



0, (38) 



from which we obtain the Green function Gq ^(^, q, k) 
ev(0^^/c^ + 9^ + fc^ + A^2 and the function w. 



Here, we have defined the 'dielectric' constant 



ev(0 = l + 



1 



(39) 



(40) 



where we go over to vortex parameters in the last equa- 
tion (m^ ^ £i) and make use of the expression ( p^ for 
the dispersive vortex elasticity £/. 

Next, we have to formulate the boundary conditions 
in terms of the zeroes of the function D^aj{k) = 0. The 
translational invariance in the y-direction allows for a 
plane- wave Ansatz axy{y) oc e"'*', while the sequence 
of 'dielectric' and 'vacuum' regions along the a;-axis, re- 
quires to match the plane waves axy{x) oc e**^^"^)"^, with 
(see (H)) 



ki = k{x < 0) = 
ko = k{0 < X < R) = k, 

kr = k{R <x)^ Jk^ - 



e 



(ev(e)-i) 



(41) 



The gauge condition Vxy ■ a^y is satisfied through the 
Ansatz 



^xy-aix)[ )e«W-)-+..) 



7 



with piecewise constant amplitudes a{x) and (3{x). The 
six coefficients are determined by the boundary condi- 
tions at X = and x — R, requiring the parallel electric 
field e|| and the magnetic field b to be continuous, see 
Appendix B. Two further conditions force the fields to 
vanish at x = ±L, L ^ oo. Requiring the determinant 
of the resulting 6x6 matrix problem to vanish, we find 



(42) 



(1 + iXO-^)ki - (e. + (Ag)-^)fc 
(l + (Ae)-2)fc, + (e, + (Ae)-2)fc 

(l + (A0-2)fc, + (e, + (A0-2)fc' 



2ikR 



For the free field, we find 



(43) 



and combining the above results for w, D, and Dq, 
Eqs. (|39|), (|4^), and (^), we can construct the function 
£)rcn|-^j-^^ g)] and obtain the formal expression for the free 
energy ( |30| ) of the vortex Casimir problem. We note that 
the cut introduced through the definitions of ki and kr in 
Eq. ( ^ ) disappears from the boundary condition D and 
hence does not contribute to the loop integral (|27|). As 
a result, the free energy expression (^^ for the Casimir 
free energy remains valid. 



C. Casimir Interaction 

In the present two dimensional situation the expression 
for the Casimir energy Eq. (|3^) per length reads 



JP'™(i?) 



I ^ r^lni^-"Ke,g)], (44) 
^'^ Jo 



and produces the Casimir force per length L 
1 aj^'™(i?) 



/ 



L OR 
-— / dq 



(45) 



{l + iXO-^)wi-{e,.. + {XO~')w 

^ il + {XO-')wr + {ey + iXO-')w _ 
(l + (A0-2)u;.-(6, + (Ae)-2)^i;^ 

We introduce the new variables 



dp / d^ 



<?\/v^ — 1 

-2N X 2 



(46) 



(l + (A0-^)g, + (6. + (A0-^). 
(l + (A0-2)s,-(6, + (A0-2)s 



Eq. ( ^ ) gives the two dimensional version for massive 
photons of the Casimir-Jorce between dielectric media, 
first derived by LifshitzcJ for the conventional 3D situa- 
tion. 

In the following we will be interested in the situa- 
tion where the 'dielectric' media are dilute (dilute vor- 
tex matter at small magnetic induction B <C <I>o/A^). 
In this case, we expand ( ^6| ) in the small correction 
tv — ^ giving the deviation of the dielectric constant 
from its vacuum value. Furthermore, if the integral is 
dominated by large frequencies ^ we can ignore the mass 
term everywhere. Under these conditions, we approxi- 
mate s ~ p and « (cv — l)(l/2p — p) + tvP, such that 
s^ — evP ~ (ev^l)(l/2p— p), while s^ + evP ~ 2p, and ob- 
tain the simplified force expression (the weak logarithmic 
dispcrion in is approximated through an appropriate 
constant) 



/ = - 



167r2 



2 foo 



dp / de 



e{2p'-i? 



c^p^ v^p^ — 1 



(47) 



In (^), ( ^ ) we have to account for three relevant 
frequency- or length scales. The first is introduced by the 
discrete nature of our superconductor: the layered struc- 
ture limits the ^-integration through the frequency cut- 
off = T^/d (this corresponds to the medium becoming 
transparent at high frequencies ^ > ^d)- The second fre- 
quency scale is introduced through the dispersion in the 
dielectric constant: at low frequencies ^ < i^a = 1/'^' the 
dielectric constant takes the form e^iS, 0) ^ 87rp/^^ 
and the dielectric response changes over to a mass renor- 
malization. The third frequency scale S,r = 1 /Re appears 
through the exponential exp(2^si?/c) > exp(2^i?£) in the 
integrand of (^) ; as ^goes beyond the exponent cuts 
off the integrand in (|4q). 

Thus, depending on the frequency ^, the media are ei- 
ther transparent, give a dielectric response, or produce a 
mass renormalization. And depending on the position of 
S.R with respect to these response regimes, we will find 
a different behavior of the Casimir force. Below, we will 
analyze the various regimes and use the pairwise summa- 
tion technique to show that the results are in agreement 
with those obtained for the van der Waals interaction 
potential in Sec. IV above. 



= v/(ev-l)+p2 + (A0-2, 
and arrive at the simplified expression 



1. Intermediate Distances X < R < d/e 

At intermediate distances X < R < d/e, we have 
Crf < Cfl such that the cutoff on the ^ integral is given 
by £,d- As the integral is dominated by large values of ^ 
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we can ignore the mass terms everywhere. The 

exponential exp(2^si?/c) restricts the p-integral to val- 
ues p < d/eR, admitting large values of p such that we 
can drop the corrections to in (^7|), e.g., p^ ~ 1 fa p^. 
Transforming p to the new variable 7 — 2R^p/c, we find 
the Casimir force on intermediate scales 



n^(ev, - 1)^ 
167r2i?2 

167r2 i?2 



d^ / d77e 



(48) 



Using the result (|l|) of the pairwise summation and in- 
serting the expression (^0|) for the dielectric constant, we 
obtain the van der Waals interaction describing the de- 
coupled limit, 



''vdW 



4eo 



T fx 



\n\-K\/d) deo \R 



(49) 



in agreement with the result (gO|) (note that the density 
p cancels in the final result for the van der Waals inter- 
action) . 



q < 1/VRK < 1/Ac. The integral dq [. . .] then 

is given by the area c/RXc times the g, ^ — > limit of the 
integrand. In the limit f — s- the mass terms are rel- 
evant and the dispersion in the dielectric constant (|40| ) 
produces the additional mass rcnormalization 



1 1 Stt 

\I- + —P^M. 



X 



A? 



(52) 



The Casimir force then decays exponentially following 
the law 



1 n-c {l-X,Mf _2fl/A. 

27r A2i? (1 + AeM)2 



(53) 



For small densities, we can expand XcM Ri 1 + AwX^p 
and obtain for the screened van der Waals interaction in 
the regime X/e < R the expression (note that the pair- 
wise summation formula (^) has to be modified for the 
exponential factor in the interaction T^vdW) 



3/2 



exp(-2i?/Ac). (54) 



2. Large Distances X,d/e < R< A/e 



At large distances, the ^-integral is cutoff on the scale 
^fl, i.e., the Casimir force is limited by the distance 
R through the exponential exp(2^si?/c) rather than by 
the transparency of the material. We still can ignore 
the mass terms in (^ as we assume that ^\ < S^r, 
providing us with a large regime for the ^-integration 
where ^ > 1/A. Transforming the energy variable ^ to 
7 = (2i?^p)/c, we obtain the following expression for the 
Casimir force in the low density limit. 



/ 



1 



87r2i?3 

19(ev - 1)2 h^c 
10247r 'W 



2 fOO f'OQ 

dj'-f'^e^'^ 
Ji 



dp 



(2p2 



1) 



16p^-\/p2 — 1 

(50) 



Following the usual scheme, this result produces the re- 
tarded van der Waals interaction 



•'vdW — 



(17l7r/256)eo T (X 



ln^(7rA/£i?) Xseo \R 



(51) 



in agreement with Eq. (|2l|). 



3. Very Large Distances X/e < R 

The situation at very large distances R > X/e = Xc 
is most conveniently analyzed in the original formulation 
( ^ ) . The mode summation is limited by the exponential 
exp(2i?x;i?), implying the following restrictions on the in- 
tegration variables ^ and q: £,/c < \/^JRXc < 1/Ac and 



VI. DISCUSSION 

We have derived the Casimir force between two bodies 
of vortex matter and have inferred from the result the 
strength of the van der Waals interaction between vor- 
tex lines via the method of pairwise summation. While 
the calculation of the Casimir force is carried out in the 
boson formulation, the van der Waals attraction is deter- 
mined in the vortex picture. The agreement between the 
results is once morejaii illustration of the equivalence of 
the two formalismsdCJ. 

The physical implications of these results lead to in- 
teresting modifications of the B-T phase diagram of lay- 
ered/anisotropic type II superconductorsB: The attrac- 
tion between the lines produces a generic vortex solid 
of density p ~ 10~^/X'^ at low temperatures. An inter- 
esting problem appearing in this context is the accurate 
determination of the entropic repulsion between the vor- 
tex lines, the latter giving an important contribution to 
the Gibbs free energy. Various approaches to this proJs,- 
lem have been discussed by Blatter and GeshkenbeinEil, 
by Volmer, Mukherji, and Nattermanncl, and by Volmer 
and SchwartzEZl. The transition from the Meissner state 
to the van der Waals vortex solid takes place via a sharp 
first-order transition at the lower critical field . The 
concurrent phase separation then can be observed in a 
decoration experiment, though very clean samples are 
required. Furthermore, the phase separation leads to the 
concentration of the vortex lines into dense regimes sep- 
arated by vortex free regions. The interaction between 
the vortex domains then is given by the Casimir force 
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calculated here, thus providing a natural reahzation of 
the physics described in this paper. 

A further remark is in place concerning the physi- 
cal character of the van der Waals interaction discussed 
here. The usual van der Waals interaction between neu- 
tral atoms arises from dipolar fluctuations in the charge 
distribution of the atoms, resulting in a force which is 
mediated through the scalar (longitudinal) potential, at 
least within the non-retarded regime at small distances. 
Comparing to the van der Waals force between the 2D 
charged bosons/vortex lines discussed here, we note that 
the elementary objects do not have an internal structure 
producing a fluctuating dipole. The van der Waals in- 
teraction then arises from current fluctuations and thus 
involves the (transverse) vector field. In principle, such 
a 'transverse' van der Waals force is also present in an 
electron gas: while the longitudinal repulsive interaction 
is screened on the Thomas-Fermi length, the transverse 
attractive interaction of the van der Waals type due to 
fluctuating currents survives at long distances (and in 
principle induces a superconducting instability). How- 
ever, the transverse van der Waals force in an electronic 
system involves the small parameter Vp/c {vp is the Fermi 
velocity), rendering the effect small. In the 2D boson 
system discussed here, the light velocity is given by the 
anisotropy of the superconductor, c = 1/e, and thus is a 
tunable parameter. 

We thank Vadim Geshkenbein, A. van Otterlo, and 
Christoph Bruder for fruitful discussions. 



The second integral I2 can be calculated in the same 
manner. 

Here, Kg denotes the 0-th order modified Bessel func- 
tions of the second kind and — + Ry. Collecting 
results and carrying out the derivatives, we obtain the 
desired result. Here, we quote the final expressions for 
the case R = (i?, 0): 



yif(i?,fc.) 



v^f{R,h) 



1 



27ra2 



27ra2 



a 
A2 



a 

a! 



Ko 



Ko 



aR 
~ 
a 
XcR 
aR 

a7 



aR 
aR\ 



oR 



(Al) 



aR 

a7 



^ ■ (A2) 



aR 
aR 
A 



The third component V^^^jR, k^) vanishes for Ry = 0. 
The important terms in (Al) and ( |A2|) are those involv- 
ing the large screening length Ac in the ar gum ent of the 
Bessel functions. The first term oc Kq in ( [A^ ) has been 
missed in Ref. H 



APPENDIX B: 2D ELECTRODYNAMICS 



APPENDIX A: PARTIAL FOURIER 
TRANSFORMS 

In this section we sketch, by way of example, the cal- 
culation of ^^.'"'(R, fcz). We start from Eq. (Ill]) and set 



a = (3 = X (here, we ignore the exponential cutoff func- 
tion due to the vortex core; the cases a = (3 = y and 
a = x, P = y follow trivially) 



1 



1 - 



(A?-A2)P 



1 + A2fc2 \^ 1 + XlK^ + \2J.2 

Performing a contour-integral over , we obtain 
T/if(R,fc.) :=/i+/2 



47rAa2 
Ac 



dfc,, 



fc2 



Y{\c,ky) 



with Y(v, w) — ^Jm^ -\- v'^w'^ and = 1 + A^fc^. Carrying 
out the integrafj over ky in Ii we find 



We sketch the main features of the 2D real time 
massive electrodynamics underlying the present work 
(we denote the mass by 1/A and restrict ourselves to 
the isotropic case). We introduce the metric g = 
diag(— 1, —1, +1) and define the gauge field a = {sLxy, at). 
In addition, we define the derivatives d" = {—Wxy,dt), 
where Wxy — {dx, dy). The field tensor /"'^ is given by 





b 



—b e. 
e 

-Cy 



y 



(Bl) 



and its dual takes the form /*" = e^^ff/s-y = 
{cy, —Cx, —b) [e"^"' is the antisymmetric tensor). The 
Lagrangian of the massive field coupled to an external 
current is given by 



C 



£ rap I a , -a 

Jap J + 2A2ff2 ° +aaj , 



(B2) 



with the coupling constant. This Lagrangian produces 
an imaginary time action 



5 [a] = Jdrd^R 



(B3) 



The transformation from reaJ-time to imaginary time is 
obtained via the formal rule^3 t —it, dt — ^ idr: o^t ^ 
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—iur- The magnetic field remains unchanged, while the 
electric field transforms to e ^ i [(V x a)^;^]^^ = ie. The 
free field real time Lagrangian (e^ — 6^) /2g'^ then goes 
over into the simple form (e^ + 6^) /2g^ = (V x a)'^/2g'^ 
within the imaginary time formalis m. 

The functional derivative of Eq. ( p32| ) provides us with 
the inhomogeneous Maxwell equations (the homogeneous 
one, 9q/*" = —9(6 + rote — 0, follows from the antisym- 
metric structure of /"^), 



(B4) 



Making use of Gauss' and Stokes' theorems for a loop 
encircling the boundary between two media, we obtain 
the boundary conditions e|| — continuous and h ~ con- 
tinuous. 

The quantization of the theory can be developed via 
the Gupta-Bleuler formalism within the Lorcntz gauge 
daO," = 0. In the case of a massive theory this method 
leads to two polarization modes, while in a massless the- 
ory the gauge invariance reduces the number of polariza- 
tion modes by one. Within the present theory, the cou- 
pling of a to the real A field in the form ia • (V x A)/ $o 
produces a gauge invariant mass term a • (a — k • (k • 
a)/fc^)/2A^(/^. As the gauge invariance is conserved in 
our theory we end up with a single polarization mode. 
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